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ABSTRACT 

We report on the study of 14 XMM-Newton observations of the magnetar SGR 1806—20 spread over a period 
of 8 years, starting in 2003 and extending to 2011. We find that in mid 2005, a year and a half after a giant 
flare (GF), the torques on the star increased to the largest value yet seen, with a long term average rate between 
2005 and 2011 of |z>| « 1.35 x 10 -11 Hz s -1 , an order of magnitude larger than its historical level measured in 
1995. The pulse morphology of the source is complex in the observations following the GF, while its pulsed- 
fraction remained constant at about 7% in all observations. Spectrally, the combination of a black-body (BB) 
and power-law (PL) components is an excellent fit to all observations. The BB and PL fluxes increased by 
a factor of 2.5 and 4, respectively, while the spectra hardened, in concordance with the 2004 major outburst 
that preceded the GF. The fluxes decayed exponentially back to quiescence with a characteristic time-scale of 
t ~ 1.5 yrs, although they did not reach a constant value until at least 3.5 years later (2009). The long-term 
timing and spectral behavior of the source point to a decoupling between the mechanisms responsible for their 
respective behavior. We argue that low level seismic activity causing small twists in the open field lines can 
explain the long lasting large torques on the star, while the spectral behavior is due to a twist imparted onto 
closed field lines after the 2004 large outburst. 


1. INTRODUCTION 

Bursting activity in magnetars (neutron stars with X-ray 
emission powered by the decay of their strong internal mag¬ 
netic fields) is usually accompanied by changes in the spectral 
and temporal properties of their persistent emission. During 
an active period, the X-ray flux of the source increases and 
can occasionally reach up to three orders of magn itude higher 
than its quiescence level (see Rea & Esposit o|2011 for a re¬ 
view). This increase is usually accompanied by spectral hard¬ 
ening (Kas pi et al.|20 03 ; Mereghet ti et al.|2015| >, changes in 
the shape of the pulse profile and the pulsed fraction (e.g., 
Muno et al.|2007 1, and variation i n the source timing proper¬ 


ties, either i n the form of a glitch (|Dib et al.|2008 | Archibald 


et al.||2013 1 or a more gradual change of the spin-down rate 


effects, however, is still not entirely clear (e.g., 

Woods et al.| 

1999, 2002"; Ng et al. 2011; Dib & Kaspi||2014 

1. These ob- 


; put tigh 

attempting to explain the X-ray emission mechanism of mag¬ 
netars. 

SGR 1806—20, one of the most prolific bursters and the last 
one to emit a Giant Flare (G F, December 2004, |Hurley et al.| 
|2005[ |Gaensler et al. |2005) >, is a perfect example of a vary¬ 
ing magnetar. Its quiescent X-ray emission shows one of the 
most erratic behaviors within the magnetar population, both 
in its timing and its spectral properties. The source power-law 
spectrum hardened gradually from an index of 2.2 in 1995 
to 1.6 in 2003 (Mereghetti et al.l 20051. The pulse profile. 


on the other hand, was consistently single pulsed during that 
time, and only changed to a double-peaked profile after the 
GF (Woods et al.|2007|. The pulse frequency and frequency 


derivative history of the source changed dramatically in the 
last twenty years (see Figure |T| of Woods et al.||20071. The 


latter decreased monotonically starting in 19997a few weeks 
after a small bursting activity episode, and reached a value 
below its 1993 historical level. It flattened out in 2002 and re¬ 


mained constant up to 2004 when strong bursting activity was 
recorded. However, a sudden and erratic change was seen 
with a hysteresis relative to that activity lasting up to a few 
months after the GF. In mid 2005, the period derivative in¬ 
creased back to the pre-2004 bursting-episode level ( |Woods| 
[et al.|2007| >. 


Here, we revisit SGR 1806—20 with emphasis on the long¬ 
term temporal and spectral variations of the source persistent 
emission. We analyze all publicly available XMM-Newton 
observations that span eight years, from 2003 to 2011. Nine 
observations took place after the December 2004 GF spread¬ 
ing over 7.5 years. The observations and data reduction are 
summarized in Section [2j Section [3] shows the evolution of 
the source temporal and spectral properties, which are then 
discussed in Section[4] 


2. OBSERVATIONS AND DATA REDUCTION 

SGR 1806—20 was observed with XMM-Newton a total of 
14 times over the span of 8 years, starting April 2003. Ten of 
these observations were performed after the December 2004 
GF. In all observations, the EPIC-pn (|S trader et al. 2001 1 > 
camera was operated in either Large Window (73-ms reso¬ 
lution) or small window (6-ms resolution) mode, using the 
thin or medium filter. All data products were obtained from 
the XMM-Newton Science Archive (XSA Qand reduced us¬ 
ing the Science Analysis System (SAS) version 13.5.0. Data 
are selected using event patterns (Ml, during only good X-ray 
events (“FLAG=0”). We applied the task c/m tp /or] to all ob¬ 
servations. This task allows for a pile-up estimate through 
the direct comparison of the fraction of the event patterns 
in a given event file to model curves from a calibration file, 
e.g., EPN_QUANTUMEF_0016.CCF for PN data. The pat¬ 
tern fraction followed the model perfectly and we concluded 


1 http://xmm.esac.esa.int/xsa/index.shtml 

2 http://xmm.esac.esa.int/sas/current/doc/epatplot/epatplot.html 
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Table 1 

Log of the XMM-Newton observations along with their timing properties 


Observation ID 

Date GTI Exposure 

z/ (error) 

PF (error)* 2 



(ks) 

(Hz) 


0148210101 b 

2003-04-03 

4.6 

0.132784(7) 

0.05 (0.02) 

0148210401'’ 

2003-10-07 

9.7 

0.1326236(4) 

0.07 (0.02) 

0205350101'’ 

2004-09-06 

43.0 

0.1323468 (7) 

0.08 (0.01) 

0164561101'’ 

2004-10-06 

17.9 

0.132325 (8) 

0.06 (0.01) 

0164561301° 

2005-03-07 

3.3 

0.13228 (2) 

0.03 (0.02) 

016456140 \ d 

2005-10-04 

31.8 

0.132156(1) 

0.07 (0.01) 

oweoosoi 11 

2006-04-04 

47.5 

0.131909(3) 

0.06 (0.02) 

040660040l e 

2006-09-10 

30.8 

0.131771 (3) 

0.06 (0.03) 

0502170301 

2007-09-26 

30.5 

0.131430 (3) 

0.07 (0.03) 

0502170401 

2008-04-02 

31.0 

0.131038(3) 

0.07 (0.02) 

0554600301 

2008-09-05 

35.0 

0.130869 (3) 

0.07 (0.03) 

0554600401 

2009-03-04 

29.6 

0.130633 (2) 

0.06 (0.02) 

0604090201 

2009-09-08 

29.0 

0.130441 (2) 

0.06 (0.03) 

0654230401 

2011-03-23 

30.0 

0.129838(1) 

0.07 (0.02) 

Notes. “ RMS pulsed fraction in 

the 1.5—10 keV energy range. Observa- 

tions also 

analyzed in: b 

Mereghetti et al. 1 2005 1 , c 

Tiengo et al. 

120051, d |Mereghetti et al. 

120071, and 

e |Esposito et al. 

(2007). 


that none were affected by pile-up. Most observations, how¬ 
ever, showed intervals of high background. In such cases, we 
excluded time periods where the background level was higher 
than 5% of the source flux. We also excluded time intervals 


assoc i ated with bursts from the source (e.g., Mereghetti et al. 
2005 1 Esposito et al. |2007| >. Finally, 
cameras from our analysis due to the 


we excluded the MOS 
poorer timing resolution 


and collecting area. Only one observation had both MOS1 
and MOS2 operating in timing mode (1.75 ms resolution, 
0164561101), 6 more observations had only MOS1 operat¬ 
ing in timing mode (0164561301, 0164561401, 0406600301, 
0406600301, 0502170301, 0502170401), while the rest were 
operated in either full-frame (2.6 s resolution, 0554600301, 
0554600401, 0654230401) or large window (0.9 s resolu¬ 
tion, 0148210101, 0148210401, 0205350101, 0604090201) 
modes. The log of the 14 XMM-Newton observations we an¬ 
alyzed is listed in Table [T] 

XMM-Newton EPIC-pn source events for all observations 
were extracted from a circle with center and radius obtained 
by running the task eregionanalyse on the cleaned event files. 
This task calculates the optimum centroid of the counts distri¬ 
bution within a given source region, and the radius of a circu¬ 
lar extraction region that maximizes the source signal to noise 
ratio. Background events are extracted from a source-free 
circle with the same radius as the source on the same CCD. 
We generated response matrix files using the SAS task rm- 
fgen, while ancillary response files were generated using the 
SAS task arfgen. The EPIC spectra were created in the energy 
range 0.5-10 keV, and grouped to have a signal to noise ratio 
of 6 with a minimum of 30 counts per bin to allow the use of 
the x 2 statistic. 


We performed our spectral analysis using XSPEC (Amaud 
1996) version 12.8.1. The photo-electric cross-sections of 


Verner et al.| ( |1996] ) and the abundances of |Wilms et al.| ( |2000l > 


are used throughout to account for absorption by neutral gas. 
All quoted uncertainties are at the ler level, unless otherwise 
noted. 



Figure 1. Frequency and frequency derivative history of SGR 1806—20 from 
mid 1993 until mid 2011. Top panel. Numbe r of bursts (per 30-da y intervals). 
The blue bars represent data collected from|Woods et al.|{2007) up to 2005 
June, and bursts reported in GCNs up to 2012. These represent bursts as seen 
with wide field-of-view instruments. We also included in red the bursts seen 
with XMM-N ewton. Mid dl e pane l. Spin frequency history. Blue dots are 
collected from |Woods et al.|j2007fr and red dots are our XMM-Newton mea¬ 
surements. The dashed ana dotted lines are fits to the frequency derivative 
from 1993 to 2000 January (z> = —1.48 x 10~ 12 Hz s -1 ), a nd 2001 January 
to 2004 April (v = -8.69 X 10" 12 Hz s _1 , Woods et al.|2007|. The dot- 
dashed line is our best fit to frequency measurements from 2005 July up to 
the last XMM-Newton observation (0 = —1.35 X 10 -11 Hz s -1 ). Bottom 
panel. Blue triangles represent instantaneous frequency derivative between 
two consecutive frequency measurements (excluding XMM-Newton data). 
Red triangles are the frequency derivative between two consecutive XMM- 
Newton observations. Note the increase in v after 2004 and the subsequent 
decrease around mid 2005. v remained more or less constant until the last 
XMM-Newton observation, with some variation around mid 2008. 


3. RESULTS 
3.1. Timing analysis 

We first applied a barycenter correction to the filtered pn 
event files. We then extracted a light curve (LC) for each of 
the 14 observations in the energy range 1.5—10 keV, at the 
128-ms and 64-ms resolution for large window and small win¬ 
dow mode observations, respectively. We ran the SAS task 
epiclccorr on these LCs to correct their count rates for back¬ 
ground and for the events lost due to various mirror inefficien- 
cie]]] We performed epoch-folding to derive an initial spin 
period for each of the observations. Then, we split each ob¬ 
servation into 5 smaller intervals, and phase-connected those 
to achieve a better spin-period measurement and smaller error. 
The results of this analysis are summarized in Table [Tj 

We used our data to extend the work of |Woods et al.f ( |2007] ). 
In that paper, the authors built a comprehensive frequency and 
frequency derivative history of SGR 1806—20 from late 1993 
up to early 2005, a few months after the GF. In Figure IT] we 
show the timing history of SGR 1806—20 extending to the last 
XMM-Newton observation (red dots), almost six and a half 
years after the GF. The blue triangles in the bottom panel of 
Figure |T| show the instantaneous frequency derivative calcu¬ 
lated between two adjacent frequency data points, excluding 
the XMM-Newton data. The red triangles show the frequency 
derivative as derived from consecutive XMM-Newton obser¬ 
vation]]] 

The combined RXTE and XMM-Newton data show some 

3 http://xmm.esac.esa.int/sas/current/doc/epiclccorr/ 

4 It is obvious here that the early XMM-Newton observations would have 
completely missed the timing noise of the source observed with RXTE in 
2004. We note, though, that the XMM-Newton data are in agreement with 
the surrounding RXTE points. 
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Figure 2. The background subtracted PPs of all XMM-Newton observations in the 1.5-10 keV, along with their different harmonic contributions. Two cycles are 
shown for clarity. The pre-GF (prior to 2005) PP has a single broad-peak morphology. It becomes more complex afterwards with contribution from at least the 
second harmonic. For clarity purposes, we omitted the error bars on the power of harmonics 3 and 4 due to their large uncertainties. 


interesting trends; the frequency derivative, almost one year 
after the GF decreased to the minimum pre-flare value (i.e., 
maximum spin-down rate; note the negative y-axis values), 
and remained at that level up to the last XMM-Newton ob¬ 
servation. The dashed and dotted lines in Figure |T] are fits 
to the frequency measurements from 1993 to 2000 January, 
v = —1.48 x 10 -12 Hz s -1 , and 2001 Jan uary to 2004 
April, v = —8.69 x 10~ 12 Hz s -1 , respectively (Woods et al. 


2007} . The dot-dashed line is our best fit to frequency mea¬ 


surements from 2005 July up to the last XMM-Newton ob¬ 
servation, v = —1.35 x 10” 11 Hz s -1 . This latter v value 
is 1.6 times larger than the one observed between 2000 and 
2004, and almost an order of magnitude larger than the histor¬ 
ical level observed prior to 2000. Deviation from the average 
after 2005 is seen around 2008 when v decreased to its lowest 
value over the course of 8 years (i.e., spin-down rate attained a 
maximum level), but recovered during the subsequent obser¬ 
vation. Unfortunately, the sparse XMM-Newton observations 
around that time did not allow for detailed modeling of this 
variation. It is interesting, however, that it happened follow¬ 
ing a modest bursting episode, as seems to be the case of the 
gradual change observed in 1999. 

For completeness, we added in the top panel of Figure [l] 
the number of bursts detected in each of the XMM-Newton 


observations (red bars) and the bursts reported in GCNs 
(Gamma-ray Coordinates Network) from 2005 up to April 
2011 (these are mostly bursts identified with the Interplane¬ 
tary Netwoi'1'0. 

Using the spin periods of the different observations, we 
epoch-folded the data to compute the pulse profiles (PPs) in 
the 1.5—10 keV energy range (Figure^. We fit the different 


PPs with a sine plus cosine function (e.g., Bildsten et al. 1997 
Younes et al. 2015[>, including contribution from four harm on - 


ics. In Figured we show the power of the different Fourier 
components relative to the total power in the signal. S imi¬ 
lar to the findings of |Gogiij et al.| ( |2002| see also |Mereghetti] 
|et al |200?) Woods et al.|2007} , we find that the PP prior to the 
December 2004 GF has a single broad-peak shape with negli¬ 
gible contribution from higher harmonics (Figure l2li. During 
and after the sixth XMM-Newton observation, i.e., 10 months 
after the GF, the PP became more complex, with contribution 
from the second harmonic, indicating a multi-peak structure, 
clearly seen in the PP (Figure EJi. 

We then derived the rms pulsed fraction (PF, e.g., Bildstenj 


5 We caution that this list of bursts,FigureflTcould be missing a few 

Itha^^ercjiotrejTorLcdm^GC^^JiowcvciM^^^jTi' 1 'd representation ol the 
I ico]uency ol the source bursting activity at a given time. 
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Figure 3. The background-subtracted PP of all XMM-Newton observations in two energy bands. Two cycles are shown for clarity. The upper rows show the 
1.5-4 keV range for each observation and the lower rows show the 4-10 keV range. See text for details. 


|et al.p997l ) of the different observations in the 1.5—10 keV 
energy range, and we estimated its error by simulating 1000 
PPs. Despite the changes in PP morphology over the 8 years, 
we find a steady PF of ~7%. This is consistent with the his¬ 
torical level measured with ASCA and BeppoSAX in 1995. 
The only change in PF was observed immediately after the 
GF, when it dropped to a minimum of 3% (observation ID 
0164561301, 2005-03-07, see also |Rea et al.|2003] >. 

Finally, we looked at the PP morphology in two different 
energy bands, 1.5—4 keV and 4—10 keV (Figure [3]). The 
two observations of September and October 2004, when the 
source flux was at a maximum, showed somewhat similar 
shapes and PFs in both bands, i.e., a single broad peak. In 
all other observations, the PPs differed between the two bands 
(the soft PP was more complex in general than the hard one), 
and the soft band PF was lower than that of the hard band. 
However, we note here that due to the low statistics of the 
data, the soft band rms PFs were not very well constrained 
and their 3cr upper limits are consistent with the hard band 
PFs. 


3.2. Spectral analysis 

We fit the pn spectra of all 14 observations simultaneously. 
We started with a simple absorbed PL model and left the PL 
indices and normalizations free to vary, while we linked the 


absorption between all spectra. We find an acceptable fit with 
a x 2 of 5402 for 5284 degrees of freedom (d.o.f.). The hy¬ 
drogen column density is JVh = (9.7 ± 0.1) x 10 22 cm -2 . 
This simple model, however, resulted in strong residuals in 
the form of a wiggle in the spectrum. We then added a BB 
component to the model (e.g., Mereghetti et akj |2005| >. This 
more complicated spectral shape provided a better lit to the 
data with a x 2 of 5180 for 5256 d.o.f. According to an F- 
test, the improvement of the PL+BB model over the single PL 
has a chance occurrence of 5.3 x 10 -22 , hence, we conclude 
that the addition of the BB is highly significant. The spectra in 
z/Fy space, are shown in Figure[4] while the spectral fit results 
are summarized in Table U 

Figure [4] shows a clear spectral evolution between the dif¬ 
ferent XMM-Newton observations, in both the PL and the BB 
components. We plotted the evolution of the different spectral 
(and temporal) parameters in Figure[5] The BB and PL fluxes 
increased by a factor of ~2.5 and 4, respectively, compared to 
their quiescent level, following the onset of the 2004 bursting 
acti vity from the source (see panel (a) and, e.g., Mereghetti] 
|et al.|[2005] >. This flux increase is clearly accompanied by 
spectral hardening in both the BB and the PL components 
(Figure [5}. 

The spectral parameters decayed back to quiescence levels 
quasi-exponentially. We fit the PL and the BB fluxes with a 
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Table 2 

BB and PL spectral parameters 


Observation ID 

kT 

r2,« 

log Fbb 

r 

log Fp L 


(keV) 

(km 2 ) 

(erg s -1 cm -2 ) 


(erg s -1 cm -2 ) 

0148210101 

0 - 1 
u '° — 0.2 

1 2 +1 ' 4 
i ' z -0.9 

-n-6±g;| 

1.3 ±0.2 

-lO.TOlg;™ 

0148210401 

0.5 ±0.2 

q 9+5.7 
°' z -2.5 

- 11.8 ± 0.6 

1 4 +0 1 
i - 4 - 0.2 

-10.63fS;“ 

0205350101 

0.90 ± 0.05 

0.7 ±0.2 

-11.23 ±0.02 

1.1 ± 0.1 

-10.41 ±0.02 

0164561101 

0.85 ± 0.06 

0.9 ±0.3 

-11.28 ±0.03 

1.2 ± 0.1 

-10.37 ±0.02 

0164561301 

0.86 ± 0.06 

1 1+ 0 - 4 

- 11.1 ± 0.2 

1 0 +0 ’ 2 
i - u -0.3 

-10.55 ±0.08 

0164561401 

0.77 ± 0.04 

1.1 ±0.3 

-11.34 ±0.02 

1.2 ± 0.1 

-10.68 ±0.03 

0406600301 

0.75 ± 0.03 

1.2 ± 0.2 

-11.28 ± 0.02 

1 0 +° 1 
lu - 0.2 

-10.83±g;gt 

0406600401 

0.75 ± 0.08 

0.5 ±0.3 

-11.29 ±0.02 

1.5 ±0.1 

-10.67 ±0.03 

0502170301 

0.74 ± 0.04 

0.9 ±0.3 

-11.48 ±0.03 

1 4 + 0 ’ 1 

-10.83 ±0.04 

0502170401 

0-69l“;“ 

1.4 ±0.3 

-11.45 ±0.02 

1.2 ± 0.2 

-10.98 ± 0.05 

0554600301 

0.65lS;“ 

1.3 ±0.3 

-11.58 ±0.03 

1.4 ±0.1 

-10.93tg : “ 

0554600401 

0.63 ± 0.04 

i O+0.4 
io -0.3 

-11.60 ±0.03 

1.4 ±0.1 

-10.99 ± 0.04 

0604090201 

0.61 ± 0.04 

1 4 + 0.5 

-11.65 ±0.03 

1.5 ±0.1 

-10.98 ± 0.04 

0654230401 

0.56 ± 0.04 

2 0+°- 7 
Z,U -0.5 

-11 64+ 0 04 

±x.uy:_o 03 

1.4 ±0.1 

-11.03 ±0.04 


Notes. “Derived by adopting an 8.7 kpc distance <|Bibby et al.|2008). Fluxes 
are in the energy range 0.5-10 keV. 



Energy (keV) 


Figure 4. vFv best fit model of all XMM-Newton observations. The sep¬ 
arate BB and PL components are shown in dashed and dotted lines, respec¬ 
tively. Notice the spectral evolution of the BB component with flux, i.e., the 
cooling trend with decreasing flux. See text for more details. 


model of the form F(t) oc exp -4 ' T (black solid lines, panels 
(f) and (h), Figure[5]l, where r is the characteristic decay time- 
scale representing a 63% decay of a given parameter back to 
quiescence. For the PL flux, we find a characteristic decay 
time-scale of r = 441 tg? 4 days. The characteristic decay 
time-scale of the BB flux was much less constrained, with 
r = 990±“i days. We find a similar decay time-scale for 
the BB temperature, kT (panel (d). Figure [5JI. A power-law 
(PL) decay trend, F(t) oc t a for the above parameters results 
in a PL decay index a = —0.65 ±0.1 for the PL flux and 
—0.4 ±0.2 for the BB flux. Finally, we derive the total energy 
emitted in each component since MJD 53126 (2004 May 01, 
i.e., the onset of the major bursting activity of 2004) assuming 
an exponential decay trend, and we find E to t, pl = 3.3 x 
10 43 ergs and E totj B b = 0.7 x 10 43 ergs for the PL and BB 
components, respectively. We find similar values when using 
the PL decay trend. 

It is worth noting that the BB area increased during the out¬ 
burst in tandem with the flux decrease from R 2 ss 0.7 km 2 


to R 2 ss 2.0 km 2 . Such a behavior has been seen before in 


a few magnetars (e.g., [Woods et al.|2004[|Israel et al.|2010|). 
Moreover, the PL to BB flux ratio reached maximum during 


the beginning of the outburst ((-Ppl/Fbb) ~ 7), where the PL 
flux increased by a factor 1.6 more than the BB flux (panel (i), 
Figure 15}. The ratio remained constant after the 2006 obser¬ 
vation ((Fpl/Fbb) ~ 4.5). 

Finally, we performed phase-resolved spectroscopy on all 
observations. We split each observation into 10 equally sep¬ 
arated bins in phase space and fit them simultaneously. We 
linked the absorbing column between the different spectra, 
and let the BB and PL spectral parameters free to vary. The 
only clear trend we observe is in the 2004 September and Oc¬ 
tober observations. We find that only the PL flux shows any 
phase variability (Mereghetti et al .||2005j l, possibly indicating 
that the PL component is responsible for most of the pulsed 
flux. In the remaining cases, unfortunately, the combination 
of worse statistics compared to the 2004 observations and the 
very low PF of the source hindered any meaningful conclu¬ 
sions. 


4. DISCUSSION 

SGR 1806—20 is one of the most active magnetars, exhibit¬ 
ing continuous bursting activity, either as strong outbursts like 
the one detected in 2004, or in sporadic groups of few iso¬ 
lated bursts (Figures [T] and [5] upper panels). SGR 1806—20 
also shows one of the most erratic behaviors in its timing 

properties (Figure [I] [Woods et al.| 2007 L _ T he very early 

observations of the source (|Kouveliotou et al.||1998[ Woods| 
et al,||2000[ Mereghetti et al.||2000 i between 1993 and 1999 
showed a steady spin-down, v = —1.48 x 10 -12 Hz s -1 . 
From early 1999 to late 2000, the spin-down rate slowly in¬ 
creased to a new steady state, v = —8.69 x 10 -12 Hz s , 
that persisted up to the outburst emitted by the source start¬ 
ing roughly in mid 2004. Thereafter, the source exhibited an 
erratic behavior leading to another dramatic change in the v 
trend ( Woods et al.|2007j . Our follow-up analysis shows that 
the spin-down of the source increased again to new long-term 
average between July 2005 and March 2011 with a new value 
of i> = —1.35 x 10 -11 Hz s -1 . 

The current SGR 1806—20 spin-down rate is an order of 
magnitude larger than its first historical value. Assuming that 
the spin-down rate is due to magnetic dipole radiation, the av¬ 
erage surface dipole magnetic field strength of a magnetar is 
B « 3.2 x 10 19 (z>/zx 3 ) 1 / 2 . where v = 1/P and P is the spin 
period of the source. An increase, therefore, of v by a fac¬ 
tor of 10 in a span of 5.5 years would imply an increase of 
3.2 in dipole magnetic field strength, which is rather unlikely. 
We explore below possible mechanisms that may have con¬ 
tributed to the torque increase on the star, hence increasing its 
spin-down rate. 

|Thompson et al~ (]2000| also see |Harding et al.|[T999| 
Thompson & Blaes 1998 [i suggested that persistent luminos- 


ity of Alfven waves and particles could generate an additional 
torque (besides the magnetic dipole one) that would affect the 
spin-down rate of a magnetar. In this model, continuous low- 
level seismic activity in a slowly spinning highly-magnetized 
neutron star will excite magnetospheric Alfven modes and in¬ 
duce a relativistic particle outflow out to a large radius. This 
particle outflow will import extra torque on the star, increasing 
its spin-down rate. According to Thompson et al.|( [2000] ), the 
spin-down torque increases by a factor of ~ (La/£mdr) 1,/2 , 
where La is the persistent luminosity of Alfven waves and 
particles, and ±mdr is the magnetic dipole luminosity (the 
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2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 



5.25 5.3 5.35 5.4 5.45 5.5 5.55 5.6 

Time (MJD, 10 4 ) 


Figure 5. Temporal and spectral histories of SGR 1806—20. Panel (a). Number of bursts in 30 da ys as seen from wide field-of-view instruments (blue bars), and 
bursts detected by XMM-Newton (red) bars. Panel (b). Frequency derivative from |Woods et al.|{2007) (blue triangles), and from XMM-Newton (red triangles). 
Panel (c). RMS PF. Panel (d). BB temperature, kT. Panel (e). BB area, R 2 . Panel (f). BB flux, Feb- Panel (g). PL index, T. Panel (h). PL flux, Fpl- Panel 
(i). PL to BB flux ratio. Panel (j). Total flux. In all panels, a dotted horizontal line represents the historical level of the given parameter measured with ASCA 
and/or BeppoSAX. The vertical dashed line is the GF epoch. Solid lines are exponential decay trends for given parameters following the 2004 outburst. See text 
for details. 
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time-average X-ray output of the SGR in quiescence derived 


using its dipole field B and spin period, Thompson & Duncan 
|l996jl. An order of magnitude increase in spin-down from 
2000 to ~mid 2005 would require an increase in the parti¬ 
cle luminosity by a factor of 100 over Lmdr (considering 
the latter to correspond to the 2000 early X-ray observations). 
An obvious effect of such a strong particle output would be 


a luminous wind nebula (Thompson et al. 

2000 

Tong et al. 

2013 see also Younes et al.j|2012). At the 8.7 
of SGR 1806—20, the angular extent of such a 

kpc distance 

wind nebula 


would be of the order of tens of arcseconds when consider¬ 
ing relativistic speeds for the particle outflow, which could 
possibly be observed with Chandra. Such a wind nebula is, 
however, not detected in SGR 1806—20 (Vigano et al.|2014j). 

A twisted magnetosphere could also impart excess torque 
on a magnetar. Magnetic stresses inside the star will cause 
motion to the footprints of the surface mag netic field lines 
caus ing these external fields to twist (Thompson et al. 2002 


Beloborodov 2009). Such a twisted configuration causes a 


slower decrease in the dipole magnetic field with increasing 
distance from the neutron star (compared to the case where no 
twist exists), resulting in a stronger B value at the light cylin¬ 
der. This excess magnetic energy is then di ssip ated through an 
increase in the spin-down rate (Thompson et Tal.|2002) i. How¬ 
ever, if the affected twisted field? are closed field lines, this 
temporal variation should also be accompanied by variation 
in the persistent X-ray spectrum, since a twisted closed field 
line wi ll have an increase in the optical depth to res onant scat¬ 
tering (Thompson et al. 2002, Beloborodov |2009] >. On the 
other hand, if the magnetic flux is dissipated through open 
field lines, it would affect the torques on the star at large radii 
without having any effect on the persistent X-ray spectrum 
(e.g., SGR J1745-2900, |Kaspi et al.|20l4] l. The spin-down 
should return to its pure dipole value as the twist angle de¬ 
creases. The time-scale for the twist energy dissipation could 
be shorter than the spin period of the source, if the twist is 
large enough that most of the energy is dissipated through 
a large magnetic reconnection event, or, more interestingly, 
for small twist angles, the twist decay time-scale could be as 
long as a few years ( Parfrey et al.||2013| |2012| |Beloborodov| 


& Thomps on|2007| l. This model could explain the continuous 
large spin-down rate of SGR 1806—20 up to mid 2011. 

We now discuss the long-term spectral evolution of the 
SGR 1806—20 X-ray emission after its 2004 major bursting 
episodcj/] SGR 1806—20, similar to many other magnetar 
sources, has a 1-10 keV spectrum well fit with the combi¬ 
nation of a thermal (BB), and a non-thermal (PL) component. 
In the magnetar model, the BB component is due to internal 
heating from the decay of the strong magnetic field, while the 
non-thermal component is due to resonant cyclotron scatter¬ 
ing of these thermal photons by the plasma in the magneto¬ 
sphere (e.g., Thompson et al. \2002) . After the 2004 bursting 
episode, the flux of the BB and PL increased, while showing 
spectral hardening (BB temperature increased and PL index, 
T, decreased). Such a fl ux-hardness relation is a prediction 
of the model presented in Thompson et al. (2002|> and a_very 


common phenomenon among magnetar sources (e.g., |Gavriil 
& Kaspf 2004] |Scholz & Kaspi|20111 >. 

One of the models usually invoked to explain the outburst 
evolution of a mag netar is of external magnetospheric origin 
(Beloborodov||20091. In this model, a twist in the magneto¬ 


sphere, if acting on a bundle of closed field lines, would in¬ 
crease their particle charge density. Returning currents from 
these field lines would hit their surface footprints, covering 
a certain area A on the surface, heating it and causing radia¬ 
tion of thermal photons. These thermal photons, in turn, are 
Compton scattered to higher energies by the plasma in the 
bundle. As the twist relaxes back to its original untwisted 
configuration, the bundle disappears gradually, and both the 
thermal and non-thermal components decrease back to their 
quiescent level. The relaxation time-scale of the bundle, 
t 10'/Lt32 < f > 1 o 1 4lii.5 S, depends primarily on the footprint 
size An .5 = (A/10 11 ' 5 ) cm 2 , the magnetic dipole moment 
M 32 = (/z/10 32 ) G cm 3 , and the electric voltage sustaining its 


plasma den sity, 4>i 0 = (<1>/10 °) V (Beloborodov 2009 


et al 


Mori 

T5 - 


. 2013). Considering a large dipole field B = 10 15 G, 
a voltage $ = 10 9 V, and the area of the BB component 
A = 47ri? 2 with R 2 ss 1 km 2 (Tablej2j>, one can roughly rec¬ 
oncile the long relaxation time for SGR 1806—20, ~ 10 8 s, 
with the one predicted by this model. However, the model 
also predicts that A decreases with decreasing flux, which we 
do not observe in the data (Figure [5}. Moreover, the PL flux 
increased by a factor of ~ 1.6 more than the BB after the on¬ 
set of the outburst (the 2004 observations), compared to the 
following observations (after October 2005) where their ratio 
was more or less constant - a behavior that does not seem to 
agree with the twisted magnetospheric model. 

Heating of the crust from re-arragement of the internal 
magnetic field, e.g., due to a sudden crack of the crust, is 
also used to reproduce the relaxation light curve of magne¬ 
tars (e.g., [Lyubarsky et aI7||2002[ |Kouveliotou et al7||2003[ 
|Pons & Rea 2012[ >. In this picture, a total energy of the or¬ 
der of ~ 10 4 ' 5 ergs is suddenly deposited into the crust, and 
re-radiated gradually in the form of thermal photons. The to¬ 
tal energy emitted from SGR 1806—20 throughout the out¬ 
burst supports these numbers. However, the decay time-scale 
for such models is expected to be much shorter (of the or¬ 
der of ~ 100 days) than the years time-scale we calculate for 
SGR 1806—20 (see discussion in |Coti Zelati et al.|2015) l. 

It is remarkable that the PF of SGR 1806—20 has remained 
at a constant value of about 7% from 1994 to 2011 (it only 
changed for a short period of time immediately after the GF, 
Rea et al.|2005|l, while all its other properties changed drasti 


cally. This could be an indication that the geometry (and loca¬ 
tion, e.g., close to the magnetic poles) of the emitting region 
has remained essentially the same for the last two decades. 

Finally, we note that there has been attempts to link the in¬ 
trinsic properties of magn etars (e.g ., v , B) with their observed 
X-ray properties (Marsde n & White|2001| Kaspi & Boydstun] 


6 The short-term evolution of SGR 1806—20 spe ctral properties after the 
2004 outburst were published injWoods et al-H2007). 


|2010[|An et al.|2012) . Such correlations assume that the mea- 

sured intrinsic properties of these sources outside of bursting 
episodes are their “true” values (assuming a dipole configu¬ 
ration). It is clear from our analysis of SGR 1806—20 (see 
also|Woods et al.|j2007j) that such assumption is not neces¬ 
sarily correct. The spin-down rate we measure for the late 
XMM-Newton observations - clearly outside of major burst¬ 
ing episodes - represents an order of magnitude variation in 
about six years. Hence, the “true” v , and by extrapolation B. 
are currently unknown for SGR 1806—20. This could also be 
true for other magnetars. The above mentioned correlations 
need to take such possibilities into account. 
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